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Про узагальнення одного методу розв'язання систем лінійних алгебраїчних рівнянь 
на випадок несиметричних матриць, що мають місце в задачах геофізики 


(Представлено академіком НАН України В. І. Старостенком) 


Будемо розглядати дійсну несиметричну матрицю 4, що діє в п-мірному евклідовому просторі К"". 
Суть методу біортогоналізації полягає в тому, що при вдалому виборі початкових векторів и, 1 у, будується 
не одна, як у методах мінімальних ітерацій К. Ланцоша чи 4А-мінімальних ітерацій, розглянутих в ||, а дві 
послідовності (и,) і |у,) взаємоспряжених векторів и, 1 у,, і - 1,2,..., п, кожна з яких знаходиться в біорто- 


гональному відношенні з однією з оболонок | (у,, А», аа (А )"у,) Р 2 Ул. | УПА А р ке дань їй Ї М 

Розгляд двох замість однієї спряжених послідовностей векторів обумовлений відомим з аналітичної гео- 
метрії фактом. Всі прості співвідношення між координатами векторів (у вигляді скалярних добутків), заданих 
в прямокутній системі, можна зберегти при переході до косокутної системи, якщо кожний вектор представити 
в двоїстому базисі двома системами взаємоспряжених координат (які в прямокутній системі між собою спів- 
падають). Задаючи вектори -- послідовні ітерації вектора несиметричною матрицею -- в двоїстому базисі 
їн.) 1 Їу,) , можемо отримати таку ж просту по структурі матрицю В, як і у випадку ітерацій симетричною 
матрицею | 1. 

Дійсно, вибудовуючи дві послідовності векторів за тричленними формулами вигляду 


ж 
Ша ЗО й о и Уа ВУ СВУ 4 Ук» УЗ п, 
ж 
ий, заці Ж 4, у зАВІМтА4 У, (1) 
і визначаючи коефіцієнти а, , 1 В, з умов біортогональності векторів и; і у,, 1-- 1,2,..., п, одержуємо 
ж 
| (4игу,) (ні,А у) 2 
М п ееривранечнім ірехаеннчни (2) 
(цу Ук) (пу Ук) 
ж 
РИ (Ай, Урі) (А узі) зд са я 
Щодо коефіцієнтів а, , 27 і краб тУччнУтчз, вони рівні між собою, що витікає з такого 


(мул) (цулз Укл) 


ланцюжка рівностей: 
(Ацз Ук) З (цу, А) 2(-а какао Пл Ай руУ зв зароГюа В рана) то 
за зара као(Ніо2 Ураі-Є какад вала КОМУ В кою (Ані зУр2) Ж 
яв ва Си) 7 (4 У) З (АГу, ма): 


Оскільки (А"у, з м, 1) З (Ай У, і) 2 (що, МУ, г) Є (цу, ЄВ камера Я ВкаУ ка) З (у, У, ), ТО, вибира- 


ючи орієнтування векторів біортогонального базису так, що (и,, У, )» 0 для будь-яких К - 1,2,..., п, маємо 
(муз Ук) 
а а Ко "К о 
Поріз Фрі о 05 Кейн В (3) 
(уз Укл) 


Так, для визначення у формулах (2) і (3) координат 0, , ,, Б,, векторів - ітерацій в біортогональному 
. « . . . ж 
базисі Ги,), (у, слід вирахувати послідовно скалярні добутки (и, |, У, |)» (ці, У,) 1 (ц,4у,)-(Аи,,у,) 
для К - 1,2,..., п (при К - Ї перший скалярний добуток не обчислюється). 

Самі вектори и, 1 У, ,, К - 1,2,...,п, знаходяться за рекурентними співвідношеннями (1). Процес об- 
числень обривається, якщо при деякому значенні К - Кк, одержуємо (и, , У, ) - 0. Вважаємо, якщо К, - п, 
то процес обчислень триває нормально, якщо (, « п -- 1, процес вироджується. Є кілька випадків виродження: 

а) и, ЗУ, 70 (двосторонній обрив); и, - 0, у, є 0 (односторонній обрив); 
6) и, 270, у, 7 0 (односторонній обрив); и, є 0, у, 20, але (и, , У, ) - 0 (тупиковий обрив). 

Односторонні і тупиковий обриви процесу обчислень являються винятковими і їх можна уникнути за до- 


помогою належного вибору початкових векторів и, і у; Двосторонній обрив аналогічний достроковому при- 


пиненню обчислень в інших методах, причини яких обговорюються в другій частині статті |. 
Виходячи з виразів (1) - - (3), можемо записати 


Лицавуцяи, Айефруцявольчни, 


(4) 





Із; : - 
Зірочкою позначена операція транспонування матриці. 


и, зо Р зо зу На й, о У вро 
для векторів з оболонки 1 (м, Ац,..., А" ц)ії 
АузВуу ну 43, 2Воу Ву КУ 
ПРИ РРО ЛОРД ОСРОНЕ (5) 
А у, аб Рв оду ЧУ,» А, ев ані кфия 
для векторів з оболонки 1. (у, А у»... (4 )" м). 


Введемо позначення: 
шов сой ТЕС Гц го ЖеЕЛеЯО, | Фо ТЕ РО ПА Мов и Те ЯК; де 








7: ПВ С КРУ 0 0 
оба ЗХ о. 30 0 0 
, 0 зб 0 0 0 
Ва Ва В 
0 0 0 жом у 2 РДА В і І 
УЮ. З збе 1 ЧО За зв 


Тоді співвідношення (4) і (5) можна записати так. АЙ «В, А У «ЦИ В. Анулюючий матрицю В 
поліном Ф , (7) будується за тричленними рекурентними формулами відповідно до схем (4) 1 (5). 


Наявність даних, які отримуються за допомогою методів, описаних в |1|, дозволяє знаходити розв'язки 
систем лінійних рівнянь 


яхебо йР ор рено (6) 

до яких зводяться багато задач геофізики. Дійсно, здійснимо перетворення 
ії ДЕ (7) 
де 0 «(ці и.,..., и, | - ортогональна матриця, складена з вектор-стовпців и,, і - 1, 2,..., п, отриманих реалі- 


зацією або методу ортогональних ітерацій, або методу мінімальних ітерацій. Помножимо рівняння (6) зліва на 
матрицю КУ" із врахуванням рівності 4(/ - (В , а також підстановки (7) отримаємо 
ВезіГЬ: (8) 
Система (8) простіша початкової (6). Особливо ця простота проявляється у випадку симетричної матриці 
А - тоді матриця В буде тридіагональною. В останньому випадку систему (8) можна розв'язати за допомо- 
гою витонченого прийому |2|, котрий можна розповсюдити на системи з несиметричною матрицею. Опису та- 
кого узагальнення немає в роботах |2-7). Він може бути важливим для застосування, тому наведемо його. 


Теорема. Якщо початкові вектори и, і У, « Ауд, Ц, У,» У; Є га структура невиродженої несиметри- 
чної матриці А: В? -» В"? 7 б. 7 бі б 
риці А: - такі, що забезпечують нормальний перебіг процесу обчислень за рекурентними 
формулами 
і 
ире4ицо еВ Мр Вища ро» Ук З 4 Ура ера Ура Ві какао» КТ 23, п, 
ж 
из зАцофуй У, З 4У Ви» 
векторів двоїстого базису Чи,), (у,) в В", де 
- (Мигзу,) 


ше 
(цу У) 


, / р ОА оте Ко ьо З 


З 
(цу 1» Ука) 
п 
то розв'язок системи лінійних рівнянь Ах - Б, який будемо шукати у вигляді Хх - Б С ій» може бути отри- 
ЗІ 
маний за п послідовних кроків наближень за схемою 
До 
1 (5, ук) кн юр (То Укл) ї ї 
х? аленбеу, ах рану, парах, хі? Й А окре НН 
(цу, у) (цу у) рак 
Доведення теореми витікає з наступних несуперечливих міркувань. Оскільки справедливі вирази 


п К-Ї 
(оте 3 бабу збо То), орі Р (худа У би, У,) є 0, ке 2 ЗМО 
ізі ізі 


1 з 


К Кк - 
то вектор х -х"" задовольняє умовам х - х" є Ци, му... М) ОХ, Ущу з У). Тому для К-ої компо- 


ненти вектора - розв'язку х в базисі (м,) можна записати 


ї р 
є из у) З (хо п) -(хох 4 Ура ЗВ раа еВ) З 


мо МД ор ЛЕ ЗРО 


- Гр» Уго 
деп з Ах- Ах? «Ь- Ах? . Звідси одержуємо, що о - Чен) ко2,3,...,п.Отже, 
(цу, У, ) 
«ет СЮ (уз Уа) 
5 -ох'абіирасічотази, (9) 
(пу У) 


що з врахуванням визначення 1-го наближення х2 відповідності з вибором початкового вектора у; - Ду 
1 0 


доводить теорему. Нев'язку вектора х"" можна вираховувати рекурентно за формулою 
С Урі) 4 


ро 
(цу Ук) 


Яна Й 
яка виводиться елементарно з посиланням на (9): 


па зро Ах? зво Ах ЗАС? ху ши пролоусео УРА кВ з ЕН 
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